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r = t — t' , even in the limit of large age t' (to < t' < t), in contrast with a 
dynamic variable in equilibrium such as the particle velocity. The equilibrium 
quantum fluctuation-dissipation theorem (QFDT) has to be modified in order to 
relate the response function Xxx(t,t') to C xx (t,t'), since this latter quantity takes 
into account even those fluctuations of the displacement which take place during 
the waiting time t w = t' — to- We describe the deviation from QFDT in terms 
of an effective inverse temperature /3 e ff.( r ;^)- The behavior of this quantity as 
a function of r for given values of T and t w is analyzed. In the classical limit it 
is shown that f3 e tt.(T,t w ) = @D(t)/[D(t) + D(t w )], where D (t) denotes the time- 
dependent diffusion coefficient. 

PACS numbers: 

05.30.-d Quantum statistical mechanics. 

05.40.+j Fluctuation phenomena, random processes, and Brownian motion. 
02.50.Ey Stochastic processes 

KEYWORDS: quantum Brownian motion, aging effects, modified quantum fluc- 
tuation-dissipation theorem. 

Corresponding author: 
Noelle POTTIER 

Groupe de Physique des Solides, Universite Paris 7, 

Tour 23, 2 place Jussieu, 75251 Paris Cedex 05, France 

Fax number: +33 1 43 54 28 78. E-mail: pottier@gps.jussieu.fr 



1 Laboratoire associe au C.N.R.S. (U.M.R. n° 7588) et aux Universites Paris 7 
et Paris 6. 

2 Laboratoire associe au C.N.R.S. (U.M.R. n° 7603) et a FUniversite Paris 6. 



2 



1. Introduction 

Out-of-equilibrium dynamics, characterized by slow relaxation and aging ef- 
fects, is currently the subject of numerous studies. Aging is by definition the 
property that the correlation function (0(t)0(t')} of an observable O at two sub- 
sequent times t' and t (£' < t) is not invariant under time translation even in 
the limit of large age t' of the system, in contrast with a dynamic variable in 
equilibrium. In other words, this correlation function depends both on the time 
difference r = t — t' and on the waiting time (or age) t w = t' . In such a non- 
equilibrium situation, the fluctuation-dissipation theorem (FDT) is violated. The 
aging phenomenon has been intensively studied, both experimentally and theoreti- 
cally, mostly in complex systems such as spin-glasses and other types of glasses [1] . 

Interestingly enough, aging can also be encountered in simpler, non disor- 
dered systems [2]- [6]. For instance, as first noticed in [2], the simplest example of 
a dynamical variable which does not reach equilibrium and exhibits aging is the 
random walk: indeed the correlation function of the displacement of a free Brow- 
nian particle at two different times depends explicitly on the two times involved, 
and not simply on their difference. In this case, the violation of the fluctuation- 
dissipation theorem can simply be characterized by a constant factor rescaling the 
temperature [2]. 

Up to now, aging effects have mostly been discussed in a classical framework 
(quantum aging in quantum spin systems have however been studied in [7]- [8]). In 
order to investigate quantum aging effects, it is interesting, to begin with, to study 
the displacement of a free quantum Brownian particle, since, as quoted above, this 
variable exhibits aging in the classical case. Free quantum Brownian motion can 
conveniently be described in the framework of the Caldeira and Leggett ohmic 
dissipation model [9]-[13]. The dissipation is introduced via a linear coupling of 
the particle to a thermal bath with a continuum of modes. In the ohmic model, 
this coupling is proportional to u> at low frequency, an assumption which allows, 
under certain conditions, to recover an equation of motion of the Langevin type. 
However, up to now, only equal time correlation functions have been computed in 
this framework and the question of aging has not yet been considered. 

The study of correlation functions at low temperature is all the more interest- 
ing since the interplay between quantum and thermal fluctuations then becomes 
effective. For instance, at low temperature, the equal time velocity correlation 
function displays a long time tail varying like — t -2 , which leads to a non-classical 
behaviour of the time-dependent diffusion coefficient D(t) as a function of time [12]. 
In the same way, non-classical effects may also be expected to take place in the 
correlation functions at two different times. This study is the aim of the present 
paper in which we compute analytically in an exact manner the two-time velocity 
and displacement correlation functions of a free quantum Brownian particle in the 
framework of the ohmic dissipation model. As a result, the two-time displacement 
correlation function displays aging whatever the temperature of the bath. The 
model allows for a discussion of the modifications of the corresponding quantum 
fluctuation-dissipation theorem (QFDT) due to aging. 
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The paper is organized as follows. In section 2, we recall briefly the results 
obtained in [2] for the two-time displacement correlation function in the viscous 
limit of the classical Langevin model and the way the classical FDT has to be mod- 
ified. This modification can be viewed as a simple rescaling of the temperature, 
the inverse temperature (3 = l/ksT being replaced by an effective inverse temper- 
ature (3 e ff, = (3/2. In section 3, remaining in the classical framework, we take the 
particle inertia into account. Then, the FDT violation may still be described by 
a rescaling of the temperature. The effective inverse temperature now effectively 
depends of the two times r and t w : /3 e ff. = Peff.( T > tw)- Interestingly enough, 
this dependence takes place via one-time quantities, namely the time-dependent 
diffusion coefficients D(t) and D(t w ): /3 e ff.( r ^™) = @D(t)/[D(t) + D(t w )]. In 
section 4, turning to the description of quantum Brownian motion, we briefly re- 
call the main features of the Caldeira and Leggett ohmic dissipation model [9]- [13]. 
Then, in section 5, we compute analytically the two-time velocity and displace- 
ment correlation functions of a free quantum Brownian particle. We analyze in 
detail the behaviour of the quantum time-dependent diffusion coefficient D®(t) 
as a function of time. This behaviour depends on the bath temperature in a way 
which allows a clear-cut distinction between a classical and a quantum regime 
for D®(t). In section 6, we discuss, in the specific case of the diffusion of a free 
quantum Brownian particle, the way the quantum fluctuation-dissipation theorem 
(QFDT) has to be modified whenever aging is present. We show how the devia- 
tions from QFDT can, for any temperature of the bath, be described in terms of an 
associated effective inverse temperature @® s {T,t w ), which can be computed from 
the quantum time-dependent diffusion coefficients D®(t) and D®(t w ). However, 
the QFDT violation cannot, except in the classical limit, be reduced to a simple 
rescaling of the bath temperature. The behavior of j3® s (r, t w ) as a function of r 
for given values of T and t w is analyzed. In Section 7, aging effects in the quantum 
diffusion problem are discussed. Finally, Section 8 contains our conclusions. 

2. Aging effects in classical Brownian motion : the viscous limit 
of the Langevin model 

In the viscous limit of the Langevin model, the one-dimensional dynamics 
of the free Brownian particle of viscosity r\ = rwy is described by the first-order 
Langevin equation 

m 1 ^ = F(t), (2.1) 

in which the Langevin force F(t) is modelized by a stationary random Gaussian 
process of zero mean and of correlation function 

(F(t)F(t')} = 2k B Tm<y 5(t - if), (2.2) 



where T is the temperature. (The symbol (...) denotes the average value over the 
realizations of the noise). 
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Let us choose some time t and consider the displacement of the particle at 
any later time t > t , as defined by 

1 /•* 

x(t)-x(t ) = / dtxF^x). (2.3) 



2.1. The two-time displacement correlation function and the dis- 
placement response function 

Let us consider the two-time displacement correlation function C xx (t,t';to) 
as defined by: 

C xx {t,t';t ) = ([x{t) - x(t )][x(t') - x{t )}). (2.4) 
It is easy to check that, for to < t' < t, 

C xx (t,t';to) = 2^(t' -t ), (2.5) 

or, in terms of the time difference r = t — t' and of the waiting time t w = t' — to, 

k T 

C xx {T,t w ) = 2 t w . (2-6) 

One recovers the fact that in this model the two-time displacement correlation 
function is proportional to the waiting time [2], [15]. 

The displacement correlation function C xx {r,t w ) is not simply a function of 
the time difference r (in this very simple model, it even does not depend on r). 
In other words, the displacement x(t) — x (to) of the free Brownian particle is not 
a stationary random process. Thus, Fourier analysis and the Wiener-Khintchine 
theorem are not applicable to this quantity. 

As for the displacement response function Xxx(t,t'), by definition, it charac- 
terizes the average displacement (x(t) —x(to)) at time t due a unit impulse of force 
taking place at a previous time t' < t. This response function is easily deduced 
from the equation of motion (2.1) in which one adds to the Langevin force a non 
random force proportional to S(t — t'). One gets 

Xxx(t,t') = Q(t-t')— , (2.7) 

m7 

where Q(t — t') denotes the unit step function. 

2.2. The modified fluctuation-dissipation theorem 

As first noticed in [2], the fluctuation-dissipation theorem (FDT) has to be 
generalized in order to handle such a case. The relation between a response func- 
tion and the associated correlation function, 

X (t,t')=(3Q(t-t')^Ml (t'<t), (2.8) 
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valid for a dynamic variable in equilibrium, is clearly not applicable to the Brown- 
ian particle displacement x(t) — x(to). As proposed in [2], one can write the actual 
relation between Xxx(t,t') and C xx (t, t'; to) with to < t' < t as 

Xxx (t, t') =(1Q(t- t') X(t, t'; t ) {tQ < t '<t), (2.9) 

where the quantity X(t, t'; to) acts as a violation factor of the fluctuation-dissipation 
theorem. Defining for the dynamic variable under study an effective temperature 
T eff .(M';*o) = l/fcsArfr. (*,*'; *o) by 

PeX.(t,t';to) = pX(t,t';to), (2.10) 
one can rewrite Eq. (2.9) as 

X**(f,0 = fa.(tJ;to) &(t - t ') dCxx ^ tQ) (to <t'< t). (2.11) 

In the viscous limit of the Langevin model, using the expressions (2.5) and 
(2.7) of C xx (t, to) and Xxx(t,t'), one gets, for any t and t' (t' < t): 

X{t,t';to) = \- (2-12) 

Thus, in this simple description of classical Brownian motion, the FDT violation 
factor is a constant [2]. Moreover, since it does not depend on T, it can simply be 
viewed as rescaling the temperature, as far as the variable of interest is concerned : 

Pea. = §• (2-13) 

3. Aging effects in classical Brownian motion : the Langevin 
model 

Actually, the full non retarded Langevin equation contains an inertial term 
and writes down 

dv s dx ,„ „ . 

m - +mlv = F(t), „=_, (3.1) 

where F(t) is the above defined Gaussian delta-correlated random force. 
3.1. The velocity correlation function 

If, at a given initial time ti the velocity takes the value Vi, at any later time 
t > U it is given by 

1 r* 

v(t)=v l e-^ t - t ^ + - dhFit^e-^-^. (3.2) 

w Ju 
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The velocity correlation function, as denned by 

C vv (t 1 ,t 2 ) = (v(ti)v(* 2 )), 

is equal to 

C vv {t u t 2 ) = (vf - ^ )e -7(t 1+ t 2 -2 fi ) + M^-il^l 
m m 

Letting ti — > — oo, the velocity v(t) reduces to its stationary part 



(3.3) 



(3.4) 



v(t) = — ( dtxFit^e 

™> J -co 



■7(*-*i) 



(3.5) 



and, correspondingly, the correlation function C vv (ti,t2) becomes the equilibrium 
correlation function 

C vv (t 1 -t 2 ) = ^e-^- t ^, (3.6) 



m 



which only depends on the time difference t\ — t 2 . Note that this latter result 
can also be obtained by taking as initial condition for the velocity a thermally 
distributed random variable Vi. 

The velocity, which equilibrates at large times, is not an aging variable. Thus, 
Fourier analysis and the Wiener-Khintchine theorem can be used equivalently in 
order to obtain the equilibrium velocity correlation function C vv (t\ — t 2 ). 

3.2. The two-time displacement correlation function and the dis- 
placement response function 

From now on, we assume that the limit ti — > — oo has been taken. The 
displacement correlation function C xx (t, t'; to) as defined by Eq. (2.4) can be ob- 
tained from the equilibrium velocity correlation function by computing the double 
integral 



C xx (t,t';to) — / dti / dt 2 C vv (t\ — t 2 ). 

J to J to 



(3.7) 



Using Eq. (3.6), one gets, for to < t' < t, 



C xx (t, t ; to) 



2(t' - 1 ) - 



1 _|_ e -7(t-*') _ e -7(*-*o) 

7 



-7(*'-*o) 



(3.8) 



or, in terms of the time difference r = t — t' and of the waiting time t w = t' — to, 



C xx (t, t w ) 



k B T 



2L 



1 + e" 



■JT 



-j(r+t w ) 



~lt^ 
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(3.9) 



Thus, as in the viscous limit, the displacement x(t) — x (to) of the free Brownian 
particle is not a stationary random process. 
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As for the displacement response function Xxx(t,t'), when one adds to the 
Langevin force in Eq. (3.1) a non random force proportional to 5(t — £'), one gets: 

Xxx (t, t') = Q(t - t') — [1 - e-T(*-*')]. (3.10) 

This quantity only depends on the time difference r = t — t'. 

3.3. The FDT violation factor and the effective temperature in the 
Langevin model 

According to Eqs. (3.8) and (3.9), the FDT violation factor in Eq. (2.9) is 
given by 

1 _ e -j(t-t') 

* « *0 = 2 _«-,(.-) - e -^-W < 3 - n > 

or, in terms of r and 

A(r,0 = — — . (3.12) 

Since X(r,t w ) does not depend on T, it can still be viewed as rescaling the tem- 
perature. The rescaling factor now depends on both times r and t w . In the limit 
of large r and t w (7T 3> 1, *yt w 3> 1), it tends towards the constant value 1/2 
which corresponds to the viscous limit of the Langevin model (Eq. (2.12)). As for 
the associated effective inverse temperature, since one has 

PeS.(r,t w )=pX(r,t w ), (3.13) 

the discussion of the behaviour of this quantity as a function of r and t w can be 
reduced to that for the FDT violation factor. 

3.4. Analysis of the violation factor and of the effective temperature 
in terms of the time-dependent diffusion coefficient 

In order to analyze at best the behaviour of X(r,t w ) and of /9 e ff.( r ^w) as a 
function of the two arguments r and t w , it is convenient to introduce the time- 
dependent diffusion coefficient D(t), which is an odd function of time defined by 

m = \j t (W)-xm 2 )- (3.i4) 

One can write 

D(t) = ^j t C xx (t,t;0). (3.15) 

For t > 0, one can use in Eq. (3.15) the above derived expression for C xx (t, t; 0). 
At any time t, one can write 

D(t) = [ duC vv {u). (3.16) 
Jo 
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In the Langevin model (3.1), using the expression (3.6) of the equilibrium 
velocity correlation function, one gets: 

7 rp 

D (t) = ^—(l- e -' rt ), t>0. (3.17) 

m7 ' 

For t > 0, D(t) increases monotonously with t. In the limit 'jt 1, it tends 
towards the usual Einstein diffusion coefficient D = hsT/trwy. 

The FDT violation factor (3.12) can be rewritten in terms of the time- 
dependent diffusion coefficients D{t) and D(t w ) as 

(3 - i8) 

This quantity is obviously smaller than 1 and, as in the viscous limit, does not 
depend on T. The associated effective inverse temperature is given by 

MrM^P ^l^ . (3.19) 

Interestingly enough, X and /? e ff. depend on r and £ w through the corresponding 
time-dependent diffusion coefficients 1 -D(r) and D{t w ). 

Actually, formula (3.18) can be understood as follows. Coming back to the 
definition (2.9) of X(t,t';to), one sees that the numerator of this quantity is pro- 
portional to the response function Xxx(t,t') while its denominator is proportional 
to the derivative dC xx (t, t'; to)/dt'. Between to and t' (i.e. during the waiting time 
t w ) : the diffusing particle does not move in the average: (x(t')) = (x(to)}. Clearly, 
in order to compute the response of the displacement to an impulse of force at 
time t' , one cannot apply the fluctuation-dissipation theorem in its standard form 
during the full time interval (to,t) (i.e. to x(t) — x(to)), but only during the re- 
stricted time interval (£', t) (i.e. to x(t) — x(t')) since in this latter case there is no 
waiting time. Indeed, doing so, that is writing 



Xxx{t,t) = /JW{t - t ) - 



(3.20) 

t =t' 



dV 

one correctly gets for the displacement response function the expression 



Xx x{t,t') = pQ{t-t') [ dhivihW)), (3.21) 

Jt' 



1 In fact, the validity of Eqs. (3.18) and (3.19) is not restricted to the simple 
(i.e. non retarded) Langevin model as described by Eq. (3.1). These formulas 
can be applied in other classical descriptions of Brownian motion in which a time- 
dependent diffusion coefficient can be defined, for instance when the motion of the 
free Brownian particle is described by a classical retarded Langevin equation. 
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or 

Xxx (t,t')=pe(t-t')D(t-t'). (3.22) 

In the Langevin model (3.1), in which the time-dependent diffusion coefficient is 
given by Eq. (3.17), the formula (3.22) for Xxx(t,t') is identical to the expression 
(3.10), as it should. As for the denominator of X(t, to), it involves the derivative 

dC xx (t,t';t ) f T . . . . 
— =/ duC vv {u), (3.23) 

that is, since C vv (u) is an even function of u (Eq. (3.6)), 

dC xx (t 7 1'; to) 



dt' 



= D(r)+D(t w ). (3.24) 



This result can be easily understood in view of the fact that fluctuations in the 
particle displacement x(t) — x(to) do take place even during the waiting time 
t w = t' — to- 

Summing up, the violation factor X(t,t';t ) and the associated effective in- 
verse temperature /3 e ff.(t, t'; to) allow to write a modified fluctuation-dissipation 
theorem relating the response function \ xx (t, t') to the derivative of a correlation 
function dC xx (t,t';to)/dt' with to < t' < t, the latter taking into account even 
those fluctuations of the particle displacement which take place during the waiting 
time. 

Note however that the parameter T e fj. = l//? e ff.) if it allows to write a modified 
FDT for the displacement of the free Brownian particle (Eq. (2.11)), cannot be 
given the full significance of a physical temperature. In particular, it does not 
control the thermalization of the free Brownian particle. Indeed, in the absence 
of potential, the thermalization, which is solely linked with the behaviour of the 
particle velocity, is effective once the limit U — > — oo has been taken. 

3.5. Discussion of the behaviour of the violation factor and of the 
effective temperature in the Langevin model 

Let us now come back to the expression (3.12) of X(r, t w ) in the classical non 
retarded Langevin model (Eq. (3.1)). Clearly, for all values of r, X is equal to 1 
when t w = and it differs from 1 when t w ^ 0. In that sense, it can be said that 
aging is always present as far as the displacement of the free Brownian particle is 
concerned. The curves representing X(r,t w ) = P e eiX T ^w)/P as a function of the 
time difference r for various non zero values of the waiting time t w are plotted 
on Fig. 1. For a given value of t w , X(r,t w ) and /3 e s.( T ^w) - like D(t) - increase 
monotonously with r towards a finite limit -^oo(^w)- 

Let us discuss the limit value of X at large values of r (7T 3> 1). In this limit, 
D{t) as given by Eq. (3.17) with t replaced by r is equal to the Einstein diffusion 
coefficient D = hsT/m^. For a given t w , the limit value of X at large values of r 
is thus 

X^) = limX(T,g = 1 . (3.25) 

7t>i a — e ' lw 
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For small values of t w ( i yt w <1), -2foo(0 does not depart very much from 1. For 
large values of t w {^t w 3> 1), X^tyj) approaches the value 1/2. The result (3.25) 
interpolates in the diffusive regime between the situation with no waiting time and 
the situation with an infinite waiting time. 

In the opposite limit, i.e. for small values of r (7T <C 1), the motion of the 
particle is ballistic, and one gets 

X ~ - 7T - . (3.26) 
For small values of t w (^t w <C 1), one has, in this regime, 

X ~ T , (3.27) 



and, for large values of 



X ~ 71-. (3.28) 



Finally, let us note that, recovering as usual the viscous limit of the Langevin 
model by letting m — > 0, 7 — > 00 with the viscosity 77 = my remaining constant, 
the factor X(r,t w ) as given by Eq. (3.12) indeed reduces to the constant value 
1/2 (Eq. (2.12)) for any r, as it should. 

4. The ohmic dissipation model 

Let us now briefly recall the main features of the ohmic dissipation model, 
which allows for a proper description of quantum Brownian motion [9]-[13]. 

4.1. The Caldeira and Leggett model 

In quantum mechanics, the equation of motion of a Brownian particle can 
be obtained by coupling it to a thermal bath constituted by a set of harmonic 
oscillators. In the Caldeira and Leggett model [11], and for a free Brownian particle 
(i.e. in the absence of an external potential, which is the situation considered here), 
the Hamiltonian of the particle-plus-bath system reads, in obvious notations, 

2 \ 2 

V v v v 

where is a real coupling constant and the last counterterm insures the trans- 
lational symmetry. From this Hamiltonian, the equations of motion for all the 
degrees of freedom can be readily written down. The bath variables once elimi- 
nated, the Heisenberg particle position operator x(t) is seen to obey the retarded 
equation of motion 

mx(t) + m I dt' k(t - t') x(t') = -m x(U) k(t - U) + F(t), (4.2) 
Jti 
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where £j denotes the initial time at which the interaction between the particle and 
the bath is switched on [11]-[12]. In Eq. (4.2), F(t) acts as a random force and k(t) 
as a retarded memory kernel. These quantities have the microscopic expressions 

k(t) = i y-=^-cosu} v t, (4.3) 
mnu u 

V 

and 

F(t) = Y J MMU)e-^ {t - U) +bl{t i )e^ t - u \ (4.4) 



4.2. The ohmic model 

In order to induce a relaxation, the oscillators of the bath must be spread 
over a continuum of frequencies. A central ingredient in the continuum limit of 
the model is the product of the density of modes of the bath, g(u>), times the 
squared coupling strength, |A(w)| 2 . One sets: 

A(u) = 27ig(u)\X(u\ 2 . (4.5) 

As far as long-time behaviours are concerned, it is sufficient to know the behaviour 
of A(u>) at low frequencies (for physical reasons, A (a;) — > when uj — > oo, as 
pictured by some cut-off function f c (u>/u> c ), where ui c is linked to the bandwidth 
of the effectively coupled bath oscillators). In the so-called ohmic dissipation 
model [11], A (a;) ~ uj for u> — > 0, and one sets 

A(u>) = m^?kof c ( — ), uj > 0, (4.6) 

UJ C 

or, equivalent ly, 

K(u) = 27/ c (-), (4.7) 

UJ C 

where K(uj) = k(t) e lujt dt denotes the Fourier transform of k(t) in the or- 
dinary sense. With a Lorentzian cut-off function f c (ui/ui c ) = uj^/(uj 2 + cu 2 ), the 
retarded memory kernel k(t) is modelized by the exponential function 

fc(t) = 7a; c e- Wc|t| . (4.8) 

In Eq. (4.8), ui~ l acts as a measure of the memory time of k{t). 

4.3. The quantum Langevin equation 

Using instead of k(t) the kernel 



k(t) = Q(t)k(t), 



(4.9) 
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the upper bound of the integral in Eq. (4.2) can be set equal to +00. This latter 
equation can then be rewritten as 



/•-t-00 

mx(t) + m dt'k(t-t')x(t') = -mx(t i )k(t-t i ) + F(t). (4.10) 
Jti 

Note that the real part of the Fourier transform K(oj) of k(t) is linked to K{uj) by 

Jfe^(w) = - K(u>). (4.11) 

Let us assume that, at time t = ti, the total density operator of the particle- 
plus-bath system is factorized in the form /9p ar t. <8pbath, with pbath = Z~ 1 e~ ,3Hbath . 
The mean value of the Langevin force F(t) is then equal to zero. As for its 
symmetrized correlation function, 

C%{t -t') = \ [(F{t)F{f) + F(t')F(t))} , (4.12) 

it depends only on the time difference t — t' and it is given by 

C%{t-t') = m — MeK(u)tw coth^e-^-t'l (4.13) 
J-00 27r 2 

The Langevin force F(t) can thus be viewed as corresponding to a stationary 
stochastic process. In the ohmic dissipation model in which K(u) is given by 
Eq. (4.7), one has, using the Lorentzian cut-off function, 

C%{t) = mfiooh r —^-^ coth ^ e"*"*. (4.14) 
J-00 2tt u 2 + u 2 c 2 v ! 

The detailed analysis of the behaviour of Cp F (t) as a function of t for any given 
temperature of the bath can be found in [12]. 

Let us now come back to the equation of motion (4.10). The r.h.s. contains, 
in addition to the Langevin force F(t), a term depending on the position x(U) of 
the particle at time t{. The kernel k(t) is supposed to have a memory time of the 
order of oj~ x and we shall study the behaviour of the system for times t such that 
Wc(t — ti) ^> 1. Then the quantity k{t — ti) is negligible, and the initial position of 
the particle becomes irrelevant in Eq. (4.10). Mathematically, this can be realized 
by letting ti — > —00. Eq. (4.10) then takes the form of the so-called quantum 
Langevin equation 



dv C +0 ° ~ dr 

m=-+m dt'k(t-t')v(t') = F(t), v=x* (4 - 15) 
^ J— 00 
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4.4. The infinitely short memory limit 

Clearly, the kernel k(t) as modelized by Eq. (4.8) admits an infinitely short 
memory limit, k(t) = 7<5(£), which is easily obtained by taking the limit u c — > oo. 
Eqs. (4.10) or (4.15) then reduce to the non retarded equation 

dv 

m— +m^v(t) = F(t), (4.16) 
at 

analogous to the classical Langevin equation (3.1), and in which memory effects 
are disregarded. Note that the corresponding limit has to be taken in the Langevin 
force correlation function (4.14), which then becomes independent of the modeliza- 
tion chosen for k(t), as it should. 

5. The two-time velocity and displacement correlation functions 
of a free quantum Brownian particle 

5.1. The two-time velocity correlation function 

The symmetrized two-time velocity correlation function is defined by 

C$ v {tiM) = \ [M*iM* 2 )+i;(*2M*i>]. 

As stated above, the Langevin force F(t) can be viewed as corresponding to a 
stationary stochastic process. Clearly, the same is true of the solution v(t) of 
Eq. (4.15), which is valid once the limit ti — > — oo has been taken. Thus, Fourier 
analysis and the Wiener-Khintchine theorem can be used to obtain the equilibrium 
velocity correlation function C^ v (ti—t 2 ), which only depends on the time difference 
ti — 1 2 . Thus, as in the classical case, the velocity of the free quantum Brownian 
particle is not an aging variable. 
As a result, one gets: 

CSL(t) = - r — J^M % u coth ^ e- (5.2) 
vvK) m J_ QO 2tt \K{uj) - iuj\ 2 2 v ; 

In the ohmic dissipation model, this expression takes the form: 



caw = 1 1 



2n (^uj c — uj 2 ) 2 + uj 2 uj 2 



hojcoth^— e~ luJt . (5.3) 



'c " 

In the case u c /"f > 4, which corresponds to a weak coupling between the particle 
and the bath 2 , one can write 



C?M — — (i - 4 -ly^ r ^ ( _L_ _ _L_) hu coth 0hL 

m v uj c ' J-^Z-k \uj 2 +uj 2 _ uj 2 + u)U 2 



(5.4) 



2 The following expressions remain valid even when uj c /"f < 4, provided that the 
appropriate analytical continuations are done. In this case, uj± has an imaginary 
part and therefore oscillations in time appear in the correlation functions. We 
shall only here consider the weak-coupling case cu c /7 > 4, since the weak-coupling 
assumption is already contained in the Caldeira and Leggett Hamiltonian (4.1). 
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where 



UJ r 



UJ± = 



1±(1-^) 1/2 

0J C 



(5.5) 



In the infinitely short memory limit u c — > oo, one simply has: 



(icu 1 



cu 2 +7" 



^coth^e"^. 



(5.6) 



At any non-zero temperature, formulas (5.4) and (5.6) can respectively be recast 
into the form 



/ / 6 U/f* 



JL/ T (f ja; _)_ JL/ T (* >a;+ ) (5.7) 



and 



<Vr ' x (/) — /•/ 



(5.8) 



where /t(*>7) stands for the series expansion 



/ T (t, 7 ) = e -7l*l +2 ^ 



-7l*l _ 



7*th. 



n|«l 
*th. 



(5.9) 



and tth. denotes a thermal time as defined by tth. = Ti/^ksT. 

The detailed analysis of the behaviour of C® v {t) as a function of t can be 
found in [12]. Let us here just quote the following important feature: there exists 
a cross-over temperature T c , which is given by T c = huj-fivkB (or T c = Try /irks in 
the infinitely short memory limit), above and below which the velocity correlation 
function exhibits two markedly different behaviours 3 . Namely, for T > T c , C® v (t) 
is positive at any time and therefore, despite the existence of quantum corrections, 
this regime may be roughly qualified as classical. For T < T c , C$ v (t) is positive 
at small times, passes through a zero and is negative at large times. This regime 
may be qualified as quantal. 

5.2. The time-dependent diffusion coefficient 

The classical definition (3.16) of the time-dependent diffusion coefficient D(t) 
can be readily extended to the present quantum model. Using the expression (5.3) 

3 This property can be proved by using the following equivalent form of fr(t, 7): 



fr(t, 7) = 7T7^th. cot 7r 7 t th . e 



-7M _ 



00 n 

2 J2 - tth ' 

n=l 



"1*1 
*th. 



With > 
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of the symmetrized velocity correlation function in the ohmic dissipation model, 
one gets for the quantum time-dependent diffusion coefficient D®(t) the expression 

D Q (t) = — — ^ ^ ft coth- — smut. 5.10 

yj m J_ 00 2tt ( 7 u c -u 2 ) 2 + uj 2 uj 2 2 y 1 

In the short-memory limit, this expression simplifies into 

D^°°(t) = 1 r ^ _L_ % coth smojt . (5.11) 

Note that formulas (5.10) and (5.11) are also valid for negative times, D®(t) and 
£)Q,w c -><x>^ being odd functions of t. 

Let us now for future purpose briefly recall the main features of the behaviour 
of D®(t) as a function of t for t > [12]. First, the limit value at infinite time 
of D®(t) is, at any non-zero temperature, the usual Einstein diffusion coefficient 
D = kBT/rwy. This result holds whatever the value, finite or not, of uj c . 

For the sake of simplicity, we restrict from now on the study of the quantum 
time-dependent diffusion coefficient to the short-memory limit u c — > oo, which 
presents the same gross features as the general case, but in which the analytical 
calculations are more simple. Formula (5.11) can (at any non-zero temperature) 
be recast into the form of the series expansion 

J*"— («) = ^ {l - «-* + ^ g \ } , t>0. (5.12) 

The quantum time-dependent diffusion coefficient thus appears as the sum of the 
classical contribution 

D u ^°°(t) = — (l-e" 7t ), t>0, (5.13) 
m 

and of a supplementary contribution due to quantum effects, which can be shown 
to be always positive. 

Above the cross-over temperature T c as defined above, £) ( 3' tJc ^°°(t) increases 
monotonously towards its limit value D = hsT/my (classical regime). Below 
the cross-over temperature (i.e. T < T c ), D^'^ ^ 00 (t) first increases, then passes 
through a maximum and finally slowly decreases towards its limit value D = 
kBT/rwy (quantum regime). Thus, in the quantum regime, i.e. when quantum 
effects are dominant, the quantum time-dependent diffusion coefficient can exceed 
its stationary value and the diffusive regime is only attained very slowly, namely 
after times t 3> t t h.- At T = 0, the quantum time-dependent diffusion coefficient 
can be expressed in terms of exponential integral functions [14] as 

D?^°°(t) = -5- [e-VMM ~ e 7t Ei(-Tf)] , t > 0. (5.14) 
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-DyL^°°(t) passes through a maximum at a time t m (T = 0) ~ 7 For 7* 3> 1, 
one has: 

D?^°°(t)~—± 7*»1- (5-15) 
7rm 7< 

The behaviour of the quantum time-dependent diffusion coefficient 
jjQ,u> c ->oo^ £ or severa i different temperatures on both sides of T c is illustrated 
on Fig. 2. Interestingly enough, these curves never intersect. One can indeed 
easily check that the slopes at the origin dD Q ' UJ ^ QO (t)/dt\ t= o = C^' Wc_>00 (t = 0) 
increase with T, as well as the asymptotic values at infinite time D = hsT/raj. 
One can also show that, at any fixed time t, the quantum time-dependent diffu- 
sion coefficient is a monotonously increasing function of the temperature 4 (Fig. 3). 
This physically expected property, which is valid in the whole temperature range 
and thus in the quantum regime (i.e. T > T c ) as well as in the classical one (i.e. 
T < T c ), will play an important role in the following. 

5.3. The two-time displacement correlation function and the re- 
sponse function 

Let us consider the symmetrized two-time displacement correlation function 
Cg(t, t'; to), as defined by 

Cg(t,t';t ) = 1 {([x(t) - x(t )][x(t') -x(t )}) + ([x(t') -x(t )][x(t) -*(*„)]>}, 

(5.16) 

where x(t) is the free quantum Brownian particle position operator. Still assuming 
that the limit U — > —00 has been taken from the beginning, C® x (t,t';to) can be 
obtained from the equilibrium symmetrized velocity correlation function C® v (t) by 
computing the double integral 



Cg(t,t';t )= I dh I dt 2 CQ v {h-t 2 ). (5.17) 

Jta J t 



However, the detailed expression of Cg x (t, t'; to) is of no use in the following anal- 
ysis, which will be seen in the next section to rely only on the expression of the 
quantum time-dependent diffusion coefficient D^(t), linked to C® x (t,t;0) by the 
relation 

D Q (t) = lj t Cg(t,t;0). (5.18) 

which is the quantum analog of the classical equation (3.15). 

As for the displacement response function Xxx(t,t'), it can be obtained from 
the equation of motion (4.15) in which one adds to the Langevin force a non 
random force proportional to 5(t — t'). One gets, in the case 5 uj c /^ > 4, and for 



4 See Section 6 (Eq. (6.12)). 

5 Formula (5.19) with u± as given by Eq. (5.5) remains valid even when 
cu c /7 < 4, provided that the appropriate analytical continuations are done. 
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any bath temperature, 



Xxx (t 1 t') = Q(t-t') x 

W17 



1 - (1 - ll)-V^ e -._(t-0 + (1 _ 47J-1/2 ^ e - w+(t - 



For well-separated time scales w c 1 and 7 1 (u c 1 7 1 ), one has simply: 



(5.19) 



Xxx (t, f ) = Q(t - t') — [1 - e -***-*'> + ^ e-^*-*')] • (5 20) 

m7 L 



In the infinitely short memory limit uj c — > 00, one recovers as expected for Xxx (t, t') 
the expression (3.10) of the Langevin model. 

6. Modified quantum fluctuation-dissipation theorem in the pres- 
ence of aging 

In the presence of aging effects, that is, when the dynamic variable of interest 
does not equilibrate with the bath, even at large times, the quantum fluctuation- 
dissipation theorem (QFDT) cannot be applied in its standard form. A modified 
form of QFDT in the non-equilibrium case allowing for the description of aging 
effects has been proposed in [7] for mean-field spin-glass models. Here, we pro- 
pose a modified form of QFDT which can conveniently be applied to the case of 
the displacement of the free quantum Brownian particle, i.e. to the problem of 
diffusion. 



6.1. Relation between the displacement response function and the 
time-dependent diffusion coefficient 

First, we shall refer to a time-dependent formulation of the equilibrium QFDT, 
namely the relation 



XBA(t,f) = —Q{t-t') 

and its reciprocal, 
dC BA (t,t') 



00 dC BA (t,t") 
dt" 



log 



tanh 



n(t" - t' 



2(3h 



(6.1) 



dt' 



= ~~ 2tt VP J dt " dt" ^- Xba ^' *") ~ XAB{t", t)] — coth 



n(t" - t') 
(3h ' 
(6.2) 

which are proved in the Appendix 6 . 

Then, let us turn to the specific problem of the diffusion of a free quantum 
Brownian particle. Since the limit t\ — > —00 has been taken from the beginning, 



6 For the sake of simplicity, the superior indexes Q in the correlation functions 
have been dropped in the general formulas. 
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the particle velocity has reached equilibrium. The equilibrium QFDT (6.1) can 
thus be applied to this dynamic variable, yielding, with obvious notations, 



2 C c 

xvx(ti,t') = — e(*i-0 j 



dt 



dt" 



-log 



tanh 



Tv(t" ~ f) 



that is 

Xvx(tl,t') 



e(*!-0 / dfc^ih-t") 



log 



tanh 



2(3% 
nit" - if) 



2(3h 



, (6.3) 



• (6-4) 



Since the displacement response function Xxx(t, t') can be expressed as the integral 



Xxx(t,t') = / Xvx(ti,t') dtx, 



(6.5) 



one gets for this quantity, making use of the identity (deduced from Eq. (5.17)) 

<)(y,M.l":l,) 



dt" 



t =t' 



duC$ v {u), 



(6.6) 



t'-t' 



the expression 



Xxx(t,t') = —Q(t-t') 



00 dt „ dcg(t,t";t ) 

— oo 



dt" 







, 7l(t"-t') 

tanh K —— ' 




-log 


t =t' 


2/3h 



(6.7) 

Note that writing formula (6.7) amounts to make use of the same argument as in 
the classical case, namely to apply QFDT as expressed by Eq. (6.1) to x(t) —x(t'). 

Since one has, in terms of the quantum time-dependent diffusion coefficient 
as defined by Eq. (5.18), 



dcg(t,t";t ) 



dt" 



t =f 



= D Q (t-t")+D Q (t" -if), 



(6.8) 



one can rewrite Eq. (6.7) as 

Xxx(t,t') = —0(t-t') J dt"[DQ(t-t")+DQ(t"-t')} 





, 7r(£ — £ ) 
tanh K —— ' 


-log 


2/3h 



-log 



that is, D Q {i) being an odd function of £, 

xxx(t,t') = —e(t-t') J dt"D^(t-t") 

One can also write equivalently 

xxx(t,t') = —e(t-t') J dt±D Q (t-t' -h) 



tanh 



iv(t" - if) 



2(3h 



(6.9) 
(6.10) 



-log 



tanh 



7Tti 



2(3h 



(6.11) 
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a formulation which clearly displays the convolution product structure of the 
displacement response function, when expressed in terms of the quantum time- 
dependent diffusion coefficient. 

Inverting this convolution product amounts to write the reciprocal rela- 
tion (6.2) for the specific problem of the diffusion of a free quantum Brownian 
particle, that is: 

= - A vp Jl ^[^(M'O-^O)]^ coth^). 

(6.12) 

Since the displacement response function Xxx does not depend on the bath temper- 
ature (Eqs. (5.19) or (5.20)), Eq. (6.12) displays the above quoted property that, 
at any fixed time t, D®(t) is a monotonously increasing function of the tempera- 
ture. In the infinitely short memory limit, taking into account the corresponding 
expression (3.10) of Xxx, one gets from Eq. (6.12): 

ft f°° 
V ; 2yrm 1 J _^ 

[S(t -t'- tje-rt-*'-*^ + e(h - t + t'^-t'-^} coth ^. (6.13) 

Eqs. (6.11) or (6.12) are the quantum generalization of the classical equation 
(3.22) and reduce to it in the classical limit, which can be easily checked by making 
use of the property 

— log | tanh ax\ ~ a ^oo ~r~ S(x). (6.14) 

Exactly like its classical analog Eq. (3.22), Eq. (6.11) allows to express the 
displacement response function Xxx(t,t') in terms of the time-dependent diffu- 
sion coefficient. However, contrary to the classical case in which Xxx (t, f) is di- 
rectly proportional to D(t — t'), in the quantum formulation Xxx(t,t') is a con- 
volution product, for the value t — t' of the argument, of the functions D®(t{) 
and — log | tanh(7r£i/2/?7i)|. This latter function is very peaked around t\ = at 
high temperature while it becomes more and more spread around this value as the 
temperature decreases. 

6.2. The modified quantum fluctuation-dissipation theorem 

Proceeding in the same way as in the classical case, let us write a modi- 
fied QFDT relating the displacement response function Xxx and the derivative 
dC^ x (t,t';to)/dt' of the correlation function C® x as defined by Eq. (5.16), with 
to < t' < £, as 

3C& (M'; to) = 
dt' 

-I ^ II *"^(*. O - X-(«*. «)] ^ coth =j£=*> , (6.15) 
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where, for short, X® stands for X®(t, to). Using the quantum analog of the 
classical equation (3.24), one gets from Eq. (6.15): 

D Q (t-t') + D Q (t' -t ) = 

"2^ VP iy ,, ^AtX)-X,At'\t) ] ^ coth^|). (6.16) 

The quantity X Q (t,t';t ) in Eqs. (6.15)-(6.16) represents the QFDT violation 
factor. 

The parameter {3® s (t,t';to) linked to X®(t,t';to) by the equation 

(3%Xt 1 t';to)=(3X < J(t,t';to), (6.17) 

which is the quantum analog of the classical equation (2.10), can be interpreted 
as an associated effective inverse temperature, defined for r = t — t' > and 
t w = t' — to > 0. In terms of this quantity, Eq. (6.16) can be rewritten as 

D Q (t-t') + D Q (t' -to) = 

~yp f dt"^[ Xxx (t, t") - Xxx (t", t)} -g— coth n{t ' ~ I } , (6.18) 

° l PeS. n PeS. h 

where for short stands for f}^ (t,t';to). 

Let us add as a comment that, like its classical analog, the parameter 
T e g = l/Pe S , if it allows to write the modified QFDT (6.18), cannot however 
be considered as a bona fide physical temperature. In particular, it does not con- 
trol the thermalization of the free quantum Brownian particle (which is effective 
here since the limit £j — > — oo has been taken from the very beginning). 

6.3. Determination of the effective temperature 

When compared with the r.h.s. of Eq. (6.12), the r.h.s. of Eq. (6.18) ap- 
pears to be formally identical to the quantum time-dependent diffusion coefficient 
D® Q (t — t') at the effective temperature T e ^ = l/ksP^ . One can thus rewrite 

Eq.° ff (6.18) as 

D Q(t _ t >) + D Q{t> _ to) = D Q {t _ (6 .i9) 



eff . 



Eq. (6.19) in turn allows for the determination of the effective temperature T e ^ 
as a function of t — t' and t' — to- Since the quantum time-dependent diffusion 
coefficient is, at any time, a monotonously increasing function of the temperature, 
Eq. (6.19) yields for (t — t', t' — t ) a uniquely defined value. 

In the classical limit in which the time-dependent diffusion coefficient is pro- 
portional to the temperature, it is easy to check that the effective inverse temper- 
ature deduced in this limit from Eq. (6.19) verifies the equation 
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in accordance with Eq. (3.19) with r = t — t' and t w = t' — t . 

Note that, except in the classical limit, f3® s as defined by Eq. (6.19) is not 
simply proportional to /3, so that the QFDT violation cannot be described by a 
simple rescaling of the bath temperature. 

6.4. Analysis of the behaviour of the effective temperature in the 
ohmic dissipation model 

In contrast with their classical counterparts, X®(t, t'; t ) and f3® s (t, t'\ t ) are 
not given by a simple ratio of time-dependent diffusion coefficients like Eq. (3.18) 
or (3.19), but have to be deduced from the implicit equation (6.19), which in 
general can only be solved numerically. Let us restrict the study, for the sake of 
simplicity, to the short-memory limit u c — > oo. Solving Eq. (6.19) amounts then 
to find the solution {3® s (t,t';t ) of the equation 

f°° cLo 1 ^ 0? s (t,t';t o )tw . , 
/ t; o o "> coth smw t — t ) = 

/ — — ~h coth^^ [smuj(t-t') + sin u(t' - t )]. (6.21) 
J.oo 2tv uj z + Y 2 

This solution is a function of the time difference r = t — t' and of the waiting time 

tw = t — to- 

The curves representing X®(r, t w ) = (3® s (r,t w )//3 as a function of r for 
various non zero values of t w are plotted on Figs. 5 to 7 for several non zero 
bath temperatures on both sides of T c . Note that, interestingly enough, Eqs. 
(6.19) and (6.21) also allow to define the effective temperature when T = 0. The 
corresponding curves are plotted on Fig. 8. 

Let us now comment upon the results obtained (Figs. 5 to 8). For all values 
of r, and at any non zero temperature, X®(t, t w ) is equal to 1 when t w = and it 
differs from 1 when t w ^ 0. Thus, like in the classical case, aging is always present 
as far as the displacement of the free quantum Brownian particle is concerned. 
For a given value of t w , X®(r,t w ) and @® s ( r ?^«;) increase monotonously with r 
towards a finite limit X^(t w ) (Figs. 5 to 7). At T = 0, (r,t w ) decreases 
monotonously with r towards zero (Fig. 8). 

Let us discuss the limit value of X® for large values of r (r 3> 7 _1 , t t h.)- For 
a given t w , one has 

Xg(t w )= lim X Q(T,t w ) = — * Q ,. r - (6.22) 

In the classical regime (i.e. T > T c ), D®(t w ) is a monotonously increasing function 
of t w , which is reflected in the fact that X^(t w ) decreases monotonously with t w 
(see Fig. 5 for T = 2T C ). In the quantum regime (i.e. T < T c ), D®(t w ) is not a 
monotonous function of t w . Consequently, the same is true of X^(t w ) (see Fig. 7 
for T = T c /2). A similar phenenomenon takes place at T = 0, in which case, for 
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any given value of r, T e ^ (r,t w ) first increases with t w , passes through a maximum 



for a value t wm 



7 , then decreases towards zero (Fig. 8). 



6.5. Discussion 

The parameter (3® s (t, t'; t ) once determined as described above, one can 
write, using Eqs. (6.18) and (6.19), 



coth 



ll(t" - t') 
fa* 



d q tQ (t-t>) = -—vp / dt"^[ Xxx (tX) -xxxV',t)] 

(6.23) 

On another hand, coming back to the expression (6.10) of Xxx(t, t'), and taking 
advantage of the fact that the displacement response function does not actually 
depend on j3 (Eqs. (5.19) or (5.20)), one may replace {3 by {3® s (t — t',t' — to) in 
this expression, which yields 



Xxx {t, t ) 



2 



/oo 
dt"D Q TQ (t - 1") 
-oo eff - 



-log 



tanh 



n{t" - if) 



2ff 



. (6.24) 



Eqs. (6.23) (resp. (6.24)) are formally similar to Eqs. (6.10) (resp. (6.12)), except 
for the fact that, due to waiting time effects, the diffusing particle is considered as 
being in contact with a bath at , an effective temperature depending on both 
time arguments t — t' and t' — to- 
Let us emphasize however that, in contrast with the sets of equations ((6.1), 
(6.2)) or ((6.10), (6.12)) which are reciprocal convolution relations, Eq. (6.24) 
cannot be viewed as the reciprocal relation of Eq. (6.23), since these two equations 
do not involve the same functions. Actually, the convolution product structure is 
lost in Eqs. (6.23) and (6.24), because of the dependence of and f3® s on 
t — t' and t' — to- As a result, while Eq. (6.23) links the displacement response 
function and D'zq i ; (t — t'), Eq. (6.24) links the displacement response 

eff ^ ^ 

function with a different function D® n , (t — t"). The dependence on the 

t% (t-t',t'-t y ' L 

running time variable t' in -D , At — t') is actually different from the t" 

6 r e « (t-t',t'-t y ' J 

dependence of D®q ^(t — t"). In particular, except for t" = t\ one has 



D Q(t - t") + D®(t" - t ) ± D% (t - t"). 



(6.25) 



7. Quantum aging and diffusion 

Keeping in mind the idea according to which aging effects in diffusion are 
basically a consequence of the fact that fluctuations in the particle displacement do 
take place even during the waiting time, a description of these effects in quantum 
diffusion can be obtained by studying the ratio 



X Q {r,t w ) = 



DQ(t) 



DQ(r) + DQ(t w )' 



r > 0, t w > 0. 



(7.1) 
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In the classical limit, this quantity is nothing but the FDT violation factor X(r, t w ) 
(Eq. (3.18)). However, outside this limit, X Q (r, t w ) is not equal to the QFDT 
violation factor X Q (r 7 t w ) = j3® s (r,t w )/(3, as defined by Eq. (6.19). 

Note that, in terms of X Q (t - t', t' - t ), Eq. (6.12) can be rewritten under 
the form of a relation between the displacement response function Xxx and the 
derivative dC^ x (t,t';t )/dt'(to < t' < t), namely 



£ *r±\X-<t, O - X^', *)] f h »th *J£=p.. (7.2) 



This equation, which generalizes the classical formula (2.9), can be viewed as an 
alternative way of writing a modified QFDT. 

The curves representing X^(r,t w ) as a function of r for various non zero 
values of t w are plotted on Figs. 9 to 11 for T = T c , T = T c /2 and T = 0. 

Let us now comment upon the results obtained (Figs. 9 to 11). For a given 
value of t w , the analysis of X®(r,t w ) closely follows the one for D®(r). Thus, in 
the classical regime (i.e. T > T c ), X^(r,t w ) increases monotonously towards its 
limit value. This is also the case when T = T c (Fig. 9). In the quantum regime 
(i.e. T < T c ), it passes through a maximum before attaining its limit (Figs. 10 
and 11). 

The limit value of X® for large values of r (r > 7 _1 ,tth.) is 

Xg(t w )= lim X<3(T,t w )= mi l nQ( , y (7-3) 

Thus, one has: 

X%(t w )=X%(t w ). (7.4) 

As indicated above, in the quantum regime, this quantity is not a monotonous 
function of t w (see Fig. 10 for T = T c /2). The limit X^(t w ) is attained after 
a time r ^> 7 _1 ,tth. (Uh. = h/2rckBT), that is, for T > T c , r ^> 7 _1 , and, for 
T < T c , t ^> tth.- This latter time becomes larger and larger as the temperature 
decreases. 

At T = 0, t t h. diverges, and a singular behaviour takes place. The diffusion 
is anomalous, as pictured by the fact that D®(t) decreases as t _1 at large time 
('jt 3> 1). For large values of r and t w (r 3> 7 _1 , t w 3> 7 _1 ), one has (Fig. 11): 

I Q (7T » 1,7*» » i) - -^r- (7-5) 

T + t w 

Note that, in contradistinction with what happens at finite temperature, the zero- 
temperature quantity possesses a singular character : as displayed by Eq. (7.5), 
the two limits r — > oo and t w — > oo do not commute. Particularly noteworthy is 
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the fact that, although in the limit £, 
Brownian particle displacement x(t) 
though diffusion is slowed down: 



-> oo, X Q tends towards 1 (Eq. (7.5)), the 
x(to) is not at thermal equilibrium, even 



([x(t) - x(t )} 2 ) ~2 log 7 (£-t ), 7(*-*o)»l- (7-6) 

nrri'j 

8. Conclusion 

This paper was devoted to the study of some aspects of a very fundamendal 
problem of non-equilibrium quantum statistical mechanics, namely the quantum 
Brownian motion of a free particle, as described in the framework of the Caldeira 
and Leggett ohmic dissipation model. We were concerned by the behaviour of 
two-time correlation functions, especially the Brownian particle velocity and dis- 
placement correlation functions C® v and C® x . Let us just here sum up some 
important features of the results. 

At a certain time ti, the particle is set in contact with a thermal bath at 
temperature T. In the limit ti — > — oo, at any finite time t the particle has been in 
contact with the bath during an infinite interval of time, and its velocity v(t) is in 
thermal equilibrium. Correspondingly, the velocity correlation function C® v (ti, £2) 
only depends on the difference t\ — £2 of the two times involved. A quantum time- 
dependent diffusion coefficient can then be defined as D®(t) = f Q C® v (u)du. For 
any given value of t, D®(t) is a monotonously increasing function of the bath 
temperature, as expected on physical grounds. 

Still assuming that the limit U — > — 00 has been taken from the very beginning, 
the situation for the particle displacement x(t) — x(to) (to < t) is completely 
different. Actually, since diffusion is going on, the variable x(t) — x(to) never 
attains equilibrium. Consequently, for any times t and t' such that to <t <t', the 
displacement correlation function C® x (t, t ) depends both on the time difference 
t = t — t' and on the waiting time or age t w = t' — t . The above analysis holds 
whatever the temperature of the bath. 

Thus, in the case of diffusion, the dynamic variable of interest, namely the 
particle displacement, does not equilibrate with the bath, even at large times, 
so that the equilibrium QFDT does not hold. In order to discuss at best how 
a modified quantum fluctuation-dissipation theorem (QFDT) can be written for 
the diffusion problem, we have chosen to refer to a time-dependent formulation of 
the equilibrium QFDT, which, in general notations, links via convolution relations 
a response function xba (t,f) and the partial derivative dCBA(t,t')/dt' of the 
associated correlation function. We have shown that it is indeed possible to write 
a modified QFDT relating the displacement response function \ xx to the partial 
derivative dC^ x (t,t';t )/dt' = D q (t) + D Q {t w ), provided that one introduces an 
associated effective inverse temperature (3^ depending on both time arguments r 
and t w . Except in the classical limit, /3® s is not simply proportional to f3 = 1/ksT, 
so that the QFDT violation cannot be reduced to a simple rescaling of the bath 
temperature. 
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The definition of heavily relies on the quantum time-dependent diffu- 
sion coefficient D® . More specifically, j3® s (r,t w ), which roughly speaking takes 
into account even those fluctuations in the particle displacement which take place 
during the waiting time, can be computed from the quantum time-dependent dif- 
fusion coefficients D®(t) and D®(t w ). Interestingly enough, the fact that {3® s is a 
uniquely defined function of r and t w finds its origin in the fact that, at any given 
time, D® s is a monotonously increasing function of the bath temperature. The 
analysis of /3® s shows that, at any bath temperature, and even at T = 0, aging is 
always present as far as the displacement of the free quantum Brownian particle is 
concerned. For any given value of t w , the effective inverse temperature /3® s (r, t w ) 
is a monotonously increasing function of r. However, the parameter = , 
albeit it allows to write a modified QFDT for the variable x(i) — x(t ) : cannot be 
given the full significance of a physical temperature, since in particular it does not 
control the thermalization of the free Brownian particle (a phenomenon which is 
effective for a particle which has been set in contact with the thermal bath for an 
infinite interval of time). 

As a natural generalization of this study, we intend, still remaining in the 
quantum framework, to consider other models for the dissipation, in which anoma- 
lous diffusive behaviours are known to take place [16], and to study the related 
two-time correlation functions and aging properties. 
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Appendix: time-dependent formulation of the equilibrium QFDT 

The equilibrium quantum fluctuation-dissipation theorem (QFDT) is usually 
written in a form which involves generalized susceptibilities and spectral densi- 
ties, which are frequency-dependent quantities. In order to characterize at best 
the modifications of QFDT due to aging effects, it is clearly more convenient 
to have at hand a formulation of the theorem in the time domain. Such a for- 
mulation establishes the link between the response function XBA(t,t') and the 
symmetrized correlation function C B A(t,t') = ±[(A(t')B(t)) + (B(t)A(t'))} or its 
derivative dCBA(t,t') /dt' . 

In order to establish this link, let us consider two quantum-mechanical ob- 
servables A and B with thermal equilibrium correlation functions verifying the 
detailed balance property 

(A(t'-ihp)B(t)) = (B(t)A(t')), (Al) 

and let us compute the contour integral 

where V is the closed contour represented on Fig. 4, with no singularities inside. 
This integral is equal to zero, according to the Cauchy theorem. In the limit 
R — > oo, one obtains, by gathering the various contributions to /, the relation 

m[(B{t)A{t')) - (A(t')B(t))] = 

/OO -t 
dt"[(A(t")B(t)) + (B(t)A(t"))\ I- , (A.3) 

~°° P smn ph — 

where the symbol vp denotes the Cauchy principal value. Taking into account the 
Kubo formula for the response function, one gets from Eq. (A3): 

9 r°° 7t i 

Wi , = _ e(t - O vp dt "C BA <t, t") ^ -^rpy- (A4) 

(3h 

Eq. ( AA) represents a formulation of QFDT in the time domain. It allows to com- 
pute the response function xba i n terms of the symmetrized correlation function 
Cba- Integrating by parts, one can recast this formula into the equivalent form 



2 r,u ^ r dC BA (t,t") 



2 r 

XBA(t,t') = —Q(t-t') J dt 



dt" 



log 



, n(t"-t') 
tanh 



2(3h 



(A5) 



In the particular case A = B, Eq. (A5) represents the quantum analog of 
the classical formula (2.8). It reduces to this latter formula in the classical limit. 
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Conversely, computing on the same contour (Fig. 4) the integral 

J = £(A(z)B(t)) jL coth dz, (A.6) 



one gets 

1 -[(A(t')B(t)) + (B(t)A(t'))] = 



? r°° ir ir(t" - t') 

"2^ VP J J f U B ® A W ~ M*"W)>] ^ coth L_. ^ (A.7) 

that is 

CBA(t,0 = -^vp dt" [ XJ m(M") ~ X AB(t",t)} jL COth ^ . (A8) 

Eq. (A. 8) represents another formulation of QFDT in the time domain. It can 
be viewed as the reciprocal of Eq. (A. 4), since it allows to compute the sym- 
metrized correlation function Cba in terms of the dissipative part of the response 
function Xba- One has also: 

dC B A(t,f) n f°° ,„ d r , „. , a ,-i i , n(t"-t') 

(-4.9) 

Like Eq. (A. 5), Eq. (A. 9) in the case A = B yields Eq. (2.8) in the classical 
limit. 

Note that, since Eqs. (A. 4) and (A. 8) are equilibrium relations, the various 
quantities into play only depend on the time differences involved, so that Eqs. (A. 4) 
and (A. 8) are in fact nothing but convolution products. In the frequency domain 
and in the case A = B, they just correspond to the usual relations 

X» = \ tanh^CV), C(u) = h coth ^ X " (a/), (A9) 

between the dissipative part x"( u ) °f the susceptibility and the Fourier transform 
C(a>) of the correlation function. 
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Figure captions 

Fig. 1 

The violation factor X(t, t w ) = /5 e ff . (t, t w ) / (3 of the classical fluctuation-dissipation 
theorem in the Langevin model plotted as a function of 7T for various values of ^yt w : 
-yt w = 0.01; -yt w = 0.1; jt w = 1; jt w = 10. 

Fig. 2 

In the short-memory limit (oj c — > 00), the quantum time-dependent diffusion coeffi- 
cient D^^c-oo^ pitted as a function of 7< for several different bath temperatures 
on both sides of T c (full lines): 
ihh. = 0.25 (T = 2T C , classical regime); 
7*th. = 0.5 (T = T c , cross-over); 
7*th. = 1 (T = T c /2, quantum regime); 
7^th. = +00 (T = 0). 

The corresponding curves for the classical diffusion coefficient D UJc ^°°(t) are plot- 
ted in dotted lines on the same figure. 

Fig. 3 

Quantum time-dependent diffusion coefficient dQi^^ 00 ^) plotted as a function 
of T/T c for several different values of jt: jt = 0.5; = 1; = 5; classical limit 
( 7 t»l). 

Fig. 4 

Integration contour for the demonstration of formulas (6.1) and (6.2) of the main 
text (or formulas (A. 5) and (A. 8) of the Appendix). 

Fig. 5 

In the short-memory limit (u c — > 00) , the violation factor X® (r, t w ) = (3® s (r, t w )//3 
of the quantum fluctuation-dissipation theorem, at bath temperature T = 2T C 
(classical regime), plotted as a function of 7T for various values of ^yt w (full lines): 
-yt w = 0.01; ^tyj = 0.1; jt w = 1; ^yt w = 10. The corresponding curves for the 
classical violation factor X(r, t w ) are plotted in dotted lines on the same figure. 

Fig. 6 

Same as above (Fig. 5), but at bath temperature T = T c (cross-over). 
Fig. 7 

Same as above (Fig. 5 or Fig. 6), but at bath temperature T = T c /2 (quantum 
regime) . 

Fig. 8 

Zero temperature case: the effective temperature (r, t w ) plotted as a function 
of 7T for various values of 7^: 7^ = 0.01; 7t TO = 0.1; ^t w = 1; ^t w = 10. 
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Fig. 9 

The ratio X Q (r,t w ) = D q (t)/[D q (t) + D Q (t w )], at bath temperature T = T c 
(cross-over), plotted as a function of 7T for various values of 7^: ^t w = 0.01; 
'jtyj = 0.1; 'jtyj = 1; 7t w = 10. 

Fig. 10 

Same as above (Fig. 9), but at bath temperature T = T c /2 (quantum regime). 
Fig. 11 

Zero temperature case: X®(t, t w ) plotted as a function of 7T for various values of 
lt w : 'jtyj = 0.01; 7t w = 0.1; ^t w = 1; ^ft w = 10; jt w = 100. 
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